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The  Response  of  a Submerged  Cylindrical  Shell 
to  an  Axially  Propagating  Acoustic  Wave 

By  G.  F.  Carrier,  Harvard  University 

1#  Introduction.  In  this  report,  we  shall  discuss  the  dynamic 
response  of  a submerged  elastic  cylindrical  shell  when  a plane 
acoustic  wave  propagates  relative  to  the  shell  as  indicated  in 
Fig.  1.  We  shall  be  concerned  only  with  the  deformation  and 
stress  in  certain  localities  where  these  quantities  are  large. 

The  results  are  found  in  a very  concise  form  and  are  of  more 
general  applicability  than  is  implied  by  the  title. 

2.  The  interaction  problem.  The  motion  of  the  shell  indicated 
in  Fig.  1 is  governed  by  the  usual  conservation  laws  as  applied 
to  the  shell  and  to  the  surrounding  fluid.  For  axially  symmetric 
deformations  of  the  shell,  the  conservation  of  momentum  require- 
ments imply  that 

EIu1^  + ohu"  + phutt  - pluj!t 

+ (Eh/R2)u  = - p(R,x,t).  (2.1) 

Here,  E,  p,  h,  R,  u are,  respectively,  the  elastic  modulus, 
density,  thickness,  radius,  and  radial  displacement  of  the  shell; 
d is  the  average  over  the  shell  thickness  of  the  axial  compres- 
sive stress  [i.e,, 

oR+h/2 

d = - i (dx/h)dr],  I = h3/12, 

; R-h/2  X ’ 

p is  the  externally  applied  over-pressure,  and  primes  denote 
differentiation  with  regard  to  the  axial  coordinate  x.  It  is  to 
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be  noted  that  we  have  retained  the  rotational  inertia  terms  in 
the  bending  theory.  Our  object,  in  doing  so,  is  to  demonstrate 
that  this  term  has  no  appreciable  effect  on  the  results  even  in 
high-speed  impact  problems.  It  should  also  be  noted  that  these 
equations  are  pertinent  to  an  unstiffened  shell.  We  shall  con- 
sider the  stiffened  shell  in  a later  section* 

The  pressure  p is  associated  with  the  acoustic  field 
surrounding  the  shell.  This  consists  of  a contribution  from  the 
incident  wave  plus  that  associated  with  the  shell  motion.  Al- 
though we  could  readily  formulate  the  problem  in  a form  which 
accounts  for  the  detailed  acoustic  field,  it  is  convenient  to 
approximate  p(R,x,t)  by 

p = p0  + pfaut(x,t)  (2.2) 

where  is  the  fluid  density,  a the  fluid  acoustic  speed,  and 
pQ  is  the  overpressure  in  the  incident  acoustic  wave.  The  Justi- 
fication of  this  approximate  formula  is  found  in  the  fact  that 
the  pressure  associated  with  the  motion  of  a plane  obstacle  nor- 
mal to  itself  is  pfaF(t),  where  F(t)  is  the  normal  velocity. 

The  equation  governing  u(x,t)  then  becomes 

EIu^  + ohuM  + phu^  " P*utt  + 

pfaUt  + = * P0S^  “ */a)  (2.3) 

where  S(t  - x/a)  is  the  conventional  step  function. 

We  shall  be  interested  in  two  basic  solutions  of  this 
equation.  The  first  is  that  for  which  u(x,t)  g w(t  - x/a). 
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This  function  must  satisfy  the  ordinary  differential  equation 
(where  5 = t - x/a) 


0^+  - )w^  + h(p  + -^)v"  + Pf.aw’ 

a4"  a 

+ (-^)w  = - p0S(O. 

R 

The  desired  solution  has  the  form 

1 P_R2  Xtf  x2^ 

-(—2 — ) [l  + a^e  + a2e  ] > 5 > 0 

W X ? Kt, 

(M!)  [a3e  3 + v 4 ] , K <0  . 

The  are  the  roots  of 

(•¥  - ^X^  + h(p  + -2-)X2  + PfaX  + = 0, 

a4-  & a2  R2 

i 

and  the  a^  are  given  by 


(2.4) 


(2.  5) 


(2.6) 


a = If  - ... . . (2.7) 

3 k/j  xj  " \ 

Ecuations  (2.5)  and  (2,7)  are  valid,  of  course,  only  if  the  Xj 
are  distinct  and  are  such  that  X^,  X2  have  negative  real  parts, 
and  X^,  have  positive  real  parts* 

An  example  of  some  interest  is  one  for  which: 


R = 8',  h = 1",  a = 5000 Vsec,  E/p  = (17,000 ’/sec)2, 
d = 25,000  #/in2,  pf/p  = 1/8,  E = 3.107#/in2* 

For  these  dimensions  (in  units,  sec"'*'),  X^  = « 650, 

Xo  = - 7500,  X . = 7300  ± 66000i.  The  coefficient,  -an , is 
c 3,t+ 


1.  These  are  chosen  so  that  u,  u',  un,  u"’,  are  continuous  at 

5=0. 
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near  unity  and  |a2!  = 0(  | X^./^  I ) > ^3!  = I 1 

p 

The  quantities  of  interest  are  w and  w"/a  since  these 
are  directly  related  to  the  hoop  and  bending  stresses,  respec- 
tively. The  orders  of  magnitude  of  these  terms  are  given  by 

'"'max  = - Po^/Eh 

lw"/a2|max  = 0[X1X2  PoR2/Eha2;l* 

It  is  now  convenient  to  turn  to  the  second  basic  solu- 
tion of  Eq.  (2.3).  We  wish,  in  fact,  to  find  the  deformation  of 
the  shell  which  is  associated  with  the  conditions  p0  = 0, 
u(x,t)  = 0 for  t < 0,  u(0,t)  = f(t),  and  u'(0,t)  = 0.  It  is 

evident  that  a combination  of  this  solution  with  the  foregoing 
will  allow  a treatment  of  practical  problems  wherein  the  shell 
is  supported  in  one  manner  or  another* 

Our  purpose  is  most  readily  accomplished  by  introducing 
the  Laplace  transform  of  u, 

— ft00  ct 

u(x,s)  = I e"  u(x,t)dt. 

J° 

The  conventional  use  of  this  transform  leads  to  the  equation 

EIu^  + (dh  - pls2)u"  + (phs2  + p^as  + Eh/R2)u  = 0. 

The  pertinent  solution  of  this  ordinary  differential  equation  is 

u(x,s)  s ----- — [qp(s)e  ^ - n-,(s)e  (2.8) 

12-11 

Here,  f(s)  is  the  transform  of  the  given  motion  u(0,t),  and 

n1>2(s)  = [M  + (M2  - n)1/2]1/2, 
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where  „ 0 

-M  = (oh  - pIs2)/2EI,  N(s)  = (phs^  + pfas  + ^g)/EI. 

R 

This  transform  would  be  rather  difficult  to  invert  with  any 
degree  of  precision  but  the  most  useful  piece  of  information  is 
readily  extracted.  The  deformation  (and  stress)  of  major  con- 
cern is  that  associated  with  the  bending  near  x = 0.  However, 


using  Eq.  (2.8),  _ i/p- 

uxx(°,s)  = -r  f 


and  this  can  be  inverted  to  give 


(2.9) 


1 /? 

uxx(°,t)  = - (ph/EI)  ' L[f (t)] , 


(2. 10) 


where  L[f(t)]  can  be  written 


2 1 -a(t—r)  I]Ik(t-T)] 

I e - 


L[f(t) ] = ft  + of  - k 


vJO 


k(t-'r) 


f (tMt 


= f^+ef+k 


f e T\l0[k(t-T)]+I^[k(t-T)])|f 1 (x)+ef (T)^dc. 


10 


(2.11) 


The  essential  feature  of  this  result  is  that  each  of  the  integrals 
in  Eqs.  (2.11)  is  positive  when  f(t)  and  f *(t)  are  positive  in 
the  time  interval  of  interest,  e.g.,  for  0 < t < tD.  Thus, 


- (EI/ph)1/2uxx(0,t) 


lies  between  the  values  f^  + ef  and  + af  where  e is  the  smaller 
root  of  N(«s)=0,  2k  is  the  positive  difference  of  the  two  roots, 
and  2a  is  the  sum  of  the  roots,  i. e. , p^a/ph.  In  the  example 
used  earlier,  a 2^  3750  sec”^  and  e ^ 650  see”^. 
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3.  Properties  of  the  foregoing  results.  Some  general  observa- 
tions associated  with  the  work  of  section  (2)  can  now  be  stated. 

We  note  first  that  the  rotational  inertia  term  has  a negligible 
effect  on  the  solution  of  the  first  problem.  In  fact,  its  in- 
clusion affects  the  size  of  and  by  about  5%  and  X-^  and  X2 
not  at  all.  In  the  latter  result,  this  term  has  no  effect  what- 
ever (rigorously).  The  deformation  at  x ^ 0 may  be  affected  by 
this  term  (must  be,  in  fact)  but  at  x = 0,  no  contribution  of 
this  item  is  present.  As  a matter  of  fact,  a simple  boundary 

1/2 

layer  type  analysis  will  demonstrate  that  only  near  x = (E/p)  t 

1/2 

and  for  x > (E/p)  t can  this  term  seriously  modify  the  result. 

We  should  also  note  that,  for  a beam  vibration  problem 
of  this  type,  the  terras  in  u and  ut  would  be  absent  and  uxx(0,t) 
is  precisely  given  by 

1/2 

uxx(°,t)  = - (ph/EI)  ut(0,t).  (3.1) 

We  should  also  note  that  the  term  representing  the  axial 

membrane  contribution  (i.e. , flu")  plays  a completely  negligible 

role.  This  is  fortunate  for  the  investigator  since,  in  many 

problems  of  interest, d will  vary  considerably  during  the  inter- 

2 

esting  time  interval  . If  d were  to  exceed  the  yield  stress,  of 
course,  this  remark  (and  the  foregoing  analysis,  as  well)  would 
be  invalid.  Note,  however,  that  during  the  deformation,  the  axial 
membrane  deformation  contributed  by  the  bending  is  such  as  to 
decrease  the  value  of  d. 

2.  As  implied  earlier,  all  of  our  remarks  apply  to  the  time 
interval  where  u and  u^.  are  positive. 
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